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ABSTRACT 

Previous papers have described the modification of the Earth Radiation Budget Satellite 
(ERBS) Attitude Determination System (ADS) to overcome the impact of onboard gyro 
de<radation and failure on attitude ground support of the mission. Two approaches were 
taken: implementing a Kalman filter in place of the batch-least-squares attitude estimator to 
account for the propagation error produced by high-noise gyro data, and modeling the ERBS 
attitude dynamics to restore rate information in the case of gyro failure. Both of these 
methods had shortcomings. In practice, the filter attitude diverged without complete sensor 
observability; and accurate dynamics modeling required knowledge of disturbance torque 
parameters that had to be determined manually. TDese difficulties have been overcome by 
improved tuning of the filter and by incorporating dynamics parameter estimation mto the 

ERBS ADS. 
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1. INTRODUCTION 


During the past 2-1/2 years, considerable effort has been devoted to overcoming the impact of 
gyro degradation and tailure in Flight Dynamics Facility ground support of the ERBS mission The ERBS 
Fine Attitude Determination System (FADS) is a batch-least-squares algorithm designed to use mo 
measurements of spacecraft motion for propagating one-orbit attitude histories, which are critical for around 
support activities such as sensor bias determination. Fine Sun sensor data for computing accurate 
single-frame attitudes are typically available for only 20 percent of one orbit. Propagation is especially 
important for yaw, which is directly observable only with the Sun sensor. At present, four of six gyro channels 
have failed completely. Prior to each failure, there was an extended period (up to 16 months) when the 
accuracy of the gyro-propagated batch attitudes was degraded by high gyro noise. Reference 1 summarizes 

the ERBS gyro performance over the mission and the impact of gyro degradation on attitude determination 
accuracy. 

Two approaches have been taken by the ERBS attitude support team to cope with the loss of accurate gyro 
data. Reference 2 describes the implementation of a Kalman filter recursive attitude estimation algorithm in 
the ERBS ADS to account for high gyro noise that caused up to 0.7-degree errors in standard batch attitudes 
The filter was found to improve pitch and roll in the case of the high gyro noise, but yaw diverged due to 
incomplete sensor observability. To restore rate information when the gyros failed completely, attitude rates 
were modeled m the ADS and used for propagation in the standard batch FADS (References 3 and 4) Results 
of the dynamics modeling were good, but several torque parameters had to be determined manually— a 
time-consuming process that made operational use of the model impractical. 

Continued efforts have improved the performance of both alternative attitude determination methods. This 
paper describes the improvements and evaluates the accuracy of each method. 


2. ERBS BACKGROUND 

The ERBS is a three-axis stabilized. Earth-oriented spacecraft, launched in 1984 into a near-circular orbit 
with an altitude of 600 km and an inclination of 57 degrees. Attitude is referenced to a geodetic coordinate 
system with pitch defined about negative orbit normal (y-axis), yaw about the local nadir vector (z-axis) and 
roll approximately along the velocity vector (x-axis). Attitude is controlled to plus or minus 1.0 degree on 
each axis. The control system used for normal flight consists of a pitch axis momentum wheel that maintains a 
strong angular momentum bias and controls pitch, two differentially driven horizon scanners mounted with 
then- axes in the y-z plane to control yaw, and electromagnetic dipoles to control roll and manage pitch axis 
angular momentum. Requirements for ground attitude determination accuracy are 0.25 degree on each axis. 
This accuracy was to be achieved using gyros, horizon scanner measurements, and digital fine Sun sensor 
data, which are usually available for only 20 percent of one 97-minute orbit. At present, the pitch and roll 
channels of both redundant gyro packages have failed. 

3. THE ERBS KALMAN FILTER 

The software modifications to convert the batch FADS to a recursive Kalman filter are minimal as many 
computations m the two algorithms are identical (Reference 2). In practice, the relative merits of each 
estimator must be considered when using them under different conditions. Batch algorithms offer robustness 
witii respect to sensor error, but are strongly dependent on the accuracy of the propagation model. A recursive 
estimator is less dependent on the accuracy of the propagation model, but is sensitive to sensor error and 
di icult to tune. Gyro-based propagation in a batch estimator gives good results as long as the gyro data are 
accurate. However, there is a point at which propagation error degrades the accuracy of batch attitudes enough 
that, even with its limitations, the Kalman filter performs better. 
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One-orbit batch and filter attitudes are compared in Figures 1 through 6, where the time span is from a period 
in the mission when the gyro data were accurate and had low noise (0.003 deg/sec root-mean-square [RMS] 
standard deviation on each axis). The reference attitude is a single frame QUatemion ESTimator (QUEST) 
solution computed using fine Sun sensor data. Single-frame QUEST attitudes computed with fine Sun sensor 
data are accurate to within 0.05 degree for roll and yaw and are considered an absolute reference for those 
axes. Due to the close alignment of the pitch axis and Sun line for the full Sun sensor coverage geometry, the 
QUEST pitch solution is based mostly on horizon scanner data and is accurate only to within about 0.2 degree. 
To simulate typical sensor observability conditions in the batch and filter runs, 80 percent of the Sun sensor 
data have been manually flagged and are not used in the estimation process. The timespan is chosen so that 
Sun sensor coverage occurs at the beg inni ng and end of the orbit. The filter pitch and roll (Figures 1 and 2) 
follow the sensor observations closely and remain within 0.3 degree of the reference QUEST pitch and roll. 
The filter yaw generally diverges when Sun sensor data are lost through the middle of the orbit. However, 
trial-and-error tuning can postpone divergence until Sun sensor data are reacquired (Figure 3). The batch 
attitude (Figures 4 through 6) diverges from the QUEST solution by up to 0.3 degree. Departures of the batch 
solution from the reference are attributed to a nonoptimal value of the epoch state, owing to the dependence on 
the less accurate horizon scanner data, and to a random walk in the gyro propagation that exists even for 
low-noise gyros. 

Figures 7 through 12 show batch and filter comparisons from a period in the mission when the pitch and roll 
gyro noise was very high (0.012 deg/sec RMS). The batch pitch and roll (Figures 7 and 8) diverge by up to 0.7 
degree. Pitch and roll propagation shows little correlation with real attitude motion. Even the batch yaw 
(Figure 9), with lower gyro noise (0.005 deg/sec RMS), is degraded through coupling of the large roll error 
into yaw. The Kalman filter (Figures 10 through 12) clearly performs better than the standard batch FADS in 
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Figure 1. One-orbit pitch comparisons for the filter solution using accurate gyro 
data for propagation (solid line) and the single frame QUEST reference 
solution (points) 

(Periods of fine Sun sensor data availability for this timespan are 
denoted by the dashed line.) 
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Figure 2. One-orbit roll comparisons for the filter solution using accurate gyro 
data for propagation (solid line) and the single frame QUEST reference 
solution (points) 



Figure 3. One-orbit yaw comparisons for the filter solution using accurate gyro 
data for propagation (solid line) and the single frame QUEST reference 
solution (points) 




Figure 4. One-orbit pitch comparisons for the batch solution using accurate gyro 
data for propagation (solid line) and the single frame QUEST reference 
solution (points) 



Figure 5. One-orbit roll comparisons for the batch solution using accurate gyro 
data for propagation (solid line) and the single frame QUEST reference 
solution (points) 
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8 6. 12. 10. 10 - 86 . 12. 10. 12-TIC = HH, hJi.SS.M 

Figure 6. One-orbit yaw comparisons for the batch solution using accurate gyro 
data for propagation (solid line) and the single frame QUEST reference 
solution (points) 



88.06.08.08 - 88 . 06 . 08. 1 0-TIC=HH. MM . SS. M 

Figure 7. One-orbit pitch comparisons for the batch solution using high noise 
gyro data for propagation (solid line) and the single frame QUEST 
reference solution (points) 

(Periods of fine Sun sensor data availability for this time span are 
denoted by the dashed line.) 
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Figure 8. One-orbit roll comparisons for the batch solution using high noise gyro 
data for propagation (solid line) and the single frame QUEST reference 
solution (points) 



Figure 9. One-orbit yaw comparisons for the batch solution using high noise gyro 
data for propagation (solid line) and the single frame QUEST reference 
solution (points) 
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88.06.OS.O8 - 88.06.08. 1 0- TI C=HH . MM. SS . M 


Figure 10. One-orbit pitch comparisons for the filter solution using high noise gyro 
data for propagation (solid line) and the single frame QUEST reference 
solution (points) 



88.06.08.08 - 88.06.08. 10-TIC=HH. MM. SS.M 

Figure 11. One-orbit roll comparisons for the filter solution using high noise gyro 
data for propagation (solid line) and the single frame QUEST reference 
solution (points) 
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Figure 12. One-orbit yaw comparisons for the filter solution using high noise gyro 
data for propagation (solid line) and the single frame QUEST reference 
solution (points) 

thi s example. However, the value of the process noise covariance matrix (the principal tuning parameter of the 
filter) that removes the divergence in yaw is inconsistent for different time spans. This difficulty may be 
overcome by determining different process noise terms for various levels of gyro noise, but it limits the 
usefulness of the filter for routine attitude determination. 

4. DYNAMICS ESTIMATION - BACKGROUND 

An alternative to gyro propagation altogether has been implemented in ERBS attitude support (References 3 
and 4). Euler’s equation for rigid-body motion was solved for the angular velocity using spacecraft control 
system telemetry data for computing the control torques and mathematical models for computing the 
disturbance torques. Using these modeled rates for propagation in the standard batch FADS , matches to within 
0.2 deg of accurate gyro-propagated solutions were obtained. However, modeled attitudes of this accuracy 
require the use of several uncertain parameters to compute the pitch axis disturbance torque. These parameters 
were identified as the spacecraft x-z product of inertia, I^, important in the pitch component of the gravity 
gradient torque; and the x and z residual dipole moments, m x and m z , important in the pitch axis magnetic 
disturbance torque. Values for these parameters were found to be inconsistent for different data spans and had 
to be determined manually for each run by trial and errror. 

To make routine use of the dynamics model practical, an automatic method of determining the uncertain 
torque parameters is necessary. The original batch FADS offers a convenient framework for estimating these 
parameters. The original FADS solves for a nine-dimension state vector consisting of the epoch attitude, gyro 
drift rate biases, and gyro scale factors. In the dynamics estimator, the three gyro scale factor state vector 
elements are replaced with the dynamics parameters m x , and m z . (Adding more state vector elements to 
the nine-dimension state would require extensive software modifications, and the gyro scale factor is not 
strictly applicable to the modeled rates anyway.) In the modified FADS, the epoch attitude estimation method 
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is unchanged; however, the epoch rate biases now correspond to the initial angular velocity in the integration 
of Euler ’s equation and are estimated slightly differently (a priori values for the three rates are taken as zero for 
roll and yaw, and the one-revolution-per-orbit rate for pitch). The only changes in the estimation are the partial 
derivatives of the attitude with respect to the new state parameters. To obtain a correct nonlinear solution of 
the dynamics parameters, the Euler solution is performed inside the FADS differential correction iteration 
loop (the new estimates of the torque parameters change the modeled rates for each iteration). 

The following is a discussion of the modifications to the FADS estimation process. The partial derivatives of 
the current attitude error, Aa, with respect to the state parameters (epoch attitude, Aa 0 , epoch angular 
velocity, Ao5 0 , and Ap comprising the dynamics parameters) are 


Aa = 
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There are two linearized differential equations governing the propagation of attitude and angular velocity 
error. 


^ = - <oAa + A(5 (2) 

^=I-'(fi'-Ol)A3 + I-'ltfAp (3) 

The primed angular momentum H' is the total angular momentum, including both that of the body 
(I(to + Aw)) and the wheel h. (The tilde denotes the antisymmetric matrix construction of a vector.) 

H' = I(w + Aw) + h (4) 

The vector AN is the error in the torque used to propagate the attitude. The columns of the torque derivative 
matrix are 


dAN 

dAIxz 


-^-(0 1- roll) 7 = (0 a - P) T 


(5) 


dAN 

dAm x 


(0 - B 3 B 2 ) t 


( 6 ) 


dAN 

dAm z 


(— B 2 B] 0) t 


(7) 


where p is the gravitational constant, r is the distance to the center of the Earth, and B is the measured 
geomagnetic field in body coordinates. 

The solutions to the attitude and angular velocity error equations have the following form 

Aa(t) = 4>(t)Aa 0 + ij)(t)Aw 0 (8) 


A6>(t) = |(t)Aw 0 + y(t)Ap 


(9) 
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If the coefficient matrices can be considered constant over a short time, the state transition matrix (<j>) is the 
matrix exponential of the coefficient matrix. 

j — 

(j) = e -(0t = cos cot [ 1 ] + (1 - cos cot) jj^jp — sin cot (10) 

The variational matrix (W) is the integral of the state transition matrix. 

t 

( sin (Dt <o (i)0) T „ 

<J)(x)dt = || w || [1] + (coscot - l)j^| - (sincot - cot)| w |p (11) 

o 

Such closed-form expressions are not now available for the angular velocity error (equation 3). These will be 
approximated by the first few terms of the power series expansion for the matrix exponential 

e At = [1] + At + A 3 t 3 + ... (12) 

where 

A = I _1 (H'-0I) ( 13 ) 

Over long time spans, when the coefficient matrices may not be constant, these equations can be solved 
recursively. The closed-form expressions and truncated series approximations above can be used over the 
short individual time steps to update the previous solution values. 

Aa ; = <J) i Aa i _ 1 + tltjAtOj (14) 


where 


AcOj = iWiAtOj.! + Y;Ap 



(15) 

(16) 


and 

t + At 

dAN; 

( 17 > 


If the partial derivative of the torque with respect to the dynamics parameters is almost constant over the short 
time step, it can be brought outside of the integral. This last expression then can be approximated as follows: 

0 0 - B 2 ‘ 

Y; « A“ 1 (e A ‘ At - [1]) I" 1 a - B 3 B, (18) 

- p B 2 0 
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The accumulated matrix for the attitude error is 


= 4>i4>i- i4»i-2 - <M>i s (19) 


The accumulated matrix for the angular velocity error is 

dAa ; aAa ; _, aAJj aAcOj aAa^, 

— — = d>: s— + — — = (t) — — 

dA( 0 o dA( 0 o dA(Oi aA( 0 o 1 dA( 0 o 


( 20 ) 


where 


aAa 0 


aA5 0 = 

(21) 

i = ?i?i- l^i-2 — 1^1 

(22) 


The accumulated matrix for the product of inertia and residual dipole moments is 

aAaj aAaj_j aAa ; aAc5j aAa^j 

— = 4>, =; — I ^ — = d> — — + xbT- 

aAp aAp aAojj aAp v aAp 1 


(23) 


Tj is computed recursively as follows: 


r i = ii r i-i + Yi 


where 


r c = [0] 


(24) 


(25) 


Combining these results, the partial derivative matrix of current attitude error with respect to the epoch 
attitude and angular velocity errors and dynamics parameters can be computed 


F i = <t>i F i-i + [ o: 


F il 


where 


F <> = {m! [ 0 ] : to] ) 


(26) 


(27) 
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5. DYNAMICS ESTIMATION - RESULTS 

To evaluate the performance of the ERBS dynamics estimation algorithm, pitch attitudes computed using the 
manual trial and error method and the dynamics estimator are compared with the gyro-propagated reference 
pitch. Figure 13 shows the modeled-to-reference pitch attitude history with the manually determined values 
of Ixz, m x , and m z . With the manual variation of parameters method, approximately 12 runs of the Euler 
solution and FADS subsystem by an experienced analyst and about 2 hours of wall clock time were required to 
achieve a match of about 0.2 degree to the reference pitch. The modeled-to-reference pitch for the dynamics 
estimation is shown in Figure 14. Automatic estimation of the dynamics parameters was accomplished in five 
differential correction iterations and in a wall clock time of only 5 minutes. Table 1 gives the values of the 
parameters for each run, together with the RMS standard deviation of the fit to reference pitch and the 
weighted observation residuals. Not only does the automatic estimation greatly increase the efficiency of the 
dynamics determination process, but it also results in a better fit to the sensor measurements than was 
accomplished manually. Roll and yaw from the dynamics estimator solution, shown in Figures 15 and 16, 
match the gyro reference to within 0.1 degree RMS. 

It is possible that the improvement in modeled pitch is due to the increased degrees of freedom in estimating 
the pitch axis disturbance torque and that the actual perturbations may be other than those modeled as gravity 
gradient and magnetic. However, these two disturbance torques are expected to be dominant for pitch. 

Knowledge of the dynamics parameters may also be used to improve the onboard control system 
performance. For example, a momentum wheel bias voltage could be uplinked to offset the constant gravity 
gradient pitch torque from the non-zero 1^ inertia product, or a dipole bias could be uplinked to offset torques 
produced by the residual magnetic dipole moment if consistent values were observed. Trend analysis of the 
dynamics parameters is also possible. Table 2 shows estimated values of I M , m x , and m z from recent mission 
data. 



Figure 13. One-orbit pitch solution using modeled rates for propagation and 

manually determined disturbance torque parameters (solid line) and 
gyro-propagated reference (points) 
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86. 12.07.04 - 86. 12.07. 06-TIC = IIH.HH.SS.tt 

Figure 14. One-orbit pitch solution using modeled rates and 

automatically determined disturbance torque parameters 
(solid line) and gyro-propagated reference (points) 
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86.12.07.04 - 86. 12.07. 06-TIC=HH . MM . SS. M 

Figure 15. One-orbit roll solution using modeled rates and automatically 
determined disturbance torque parameters (solid line) and 
gyro-propagated reference (points) 
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86.12.07.84 - 86. 12.07. e6-TIC=IIH.HM.SS.M 


Figure 16. One-orbit yaw solution using modeled rates and automatically 
determined disturbance torque parameters (solid line) and 
gyro-propagated reference (points) 


Table 1 . Dynamics Estimation Comparisons for 861 207 


METHOD 

l« (kg— m2) 

m x (ATm2) 

m z (ATm2) 

FIT TO 
REFERENCE 
PITCH (RMS) 

WEIGHTED 
OBSERVATION 
RESIDUALS (RMS) 

MANUAL 

-20.00 

3.50 

0.50 

0.12 

0.901 

DYNAMICS 

ESTIMATOR 

-20.47 

3.02 

0.26 

0.06 

0.816 

GYRO 

REFERENCE 

— 

— 



0.776 


Table 2. Dynamics Parameter Trends 


DATE 

I xz (kg-m2) 

m x (ATm2) 

m z (ATm2) 

RESIDUALS (RMS) 

920311 

-42.6 

1.71 

-2.26 

0.91 

920320 

-48.5 

-0.51 

1.03 

1.21 

920325 

-52.1 

0.42 

2.05 

1.35 

920401 

-42.8 

0.93 

2.64 

1.13 


6. CONCLUSIONS 

Two strategies to overcome the loss of accurate gyro data in ground attitude determination support of the 
ERBS mission have been implemented. Both methods involve a modification of the existing ADS, taking 
advantage of existing software to minimize development effort. 
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Although the Kalman filter performs better than the batch estimator for very high gyro noise, the filter yaw 
solution is very sensitive to the process noise. Divergence in the filter yaw may be reduced by appropriate 
tuning of the process noise terms, but the tuning was found to be inconsistent even for time spans with the 
same gyro noise levels. This inconsistency in yaw behavior makes routine use of the filter impractical. An 
adaptive Kalman filter that determines the process noise automatically would be more suitable for the 
application to ERBS attitude determination, where the propagation noise has varied widely over the course of 
the mission. 

The dynamics estimator is a viable solution to the problem of gyro failure for ERBS attitude ground support. 
Modeled rates can be used to propagate one-orbit attitudes to an accuracy within the 0.25-degree requirement 
with no a priori knowledge of disturbance torque parameters. Values of the dynamics parameters determined 
in the estimation process are also useful for analysis of control system performance. The dynamics estimator 
is currently being evaluated for operational use to restore full Flight Dynamics Facility attitude ground 
support of the ERBS mission. 
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